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Electricity consumption is increasing gradually and this trend will continue 
in the future. In addition, rapid network control systems using the resources 
offered by power electronics and control microelectronics have been 
recently studied and developed, and are currently in normal application for 
some, for others, in pilot applications or as prototypes. This paper attempts 
to show that these systems are referred to by the general acronym flexible 
alternative current transmission systems (FACTS) similarly dethroned the 
traditional systems while offering better solutions and solving the energy 
quality problem such as the hybrid system (unified power flow controller 
(UPFC), or universal phase shifter regulator (UPSR)) which opens up new 
perspectives for more efficient operation of networks by continuous and 
rapid action on the various parameters of the network (voltage, phase shift, 
and impedance); thus, the power transits will be better controlled and the 
voltages better held, which will make it possible to increase the stability 
margins or tend towards the thermal limits of the lines. In this work, we used 
a classic control (PI-decoupled) and others while offering more flexibility of 
control thanks to the development of strategies identification/control based 
on generalized predictive control (GPC) with neural network to ensure 
robust control with advanced algorithms. 
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1. INTRODUCTION 


As the demand for power increases, the operation of the power system to operate within power 
quality standards becomes more complex and therefore less secure. The advent of flexible alternating current 
transmission system (FACTS) technologies has coincided with a major restructuring of power generation; 
they can provide significant benefits by increasing the transmission capacity of the system as well as the 
flexibility and speed of power flow control. There are several types of FACTS, and the choice of the 
appropriate device depends on the goal to be achieved [1]-[3]. The unified power flow controller is one of the 
most effective FACTS equipment systems for improving power system safety; however, the performance of 
the UPFC is highly dependent on the parameter settings of the equipment in the system. The utilization of 
static VAR compensator (SVC) converters in power flow control for AC electric power transmission systems 
is an established practice, the UPFC is a union of two SVCs one is a static synchronous compensator 
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(STATCOM) shunt device and the other is a static synchronous series compensator (SSSC) series 
component, the two devices can operate independently, they act on the control of the voltage and/or the 
impedance of the line by injecting the amount of active and/or reactive power required, or simultaneously. 
The unified power flow controller (UPFC) combines both series compensation and parallel compensation, it 
is the most powerful device in power flow control, where it allows for adjustment of the three network 
quantities namely, voltage, angle load, and line impedance. The basic structure of a UPFC connected to the 
electrical power network is represented by Figure | [4], [5]. It consists of two three-phase inverters, one 
connected in series with the network via three single-phase transformers whose primaries are connected, in a 
star, and the other connected in parallel with the network via a three-phase transformer. The two inverters are 
interconnected by a DC bus represented by capacitor C. 


val Inverter I Inverter II 


Eig 


Figure 1. The basic scheme of a UPFC 


2. MODELING OF THE UPFC 

The two-level UPFC is characterized by three systems of equations which are detailed in the 
following paragraphs [6]: i) the dynamic equations of the series compensator, ii) the dynamic equations of 
the parallel compensator, and iii) the equations of the continuous circuit. The application of Kirchhoff's laws 
to meshes of the circuit of Figure 2 and the transformation of PARK, the current in the transmission line can 
be described by (1) to (3). 


digg _ ä r, 1 
u ~ Lis Iisa tT (Vsa — Vea — Vra ) 
disg r, 1 (1) 
dt Ulsa — L lsq + 7 Usa = Veq T Urq) 
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Figure 2. Equivalent circuit of UPFC system 


3. CONTROLLER DESIGN 

This study provides new insights into what was designed to improve the performance of the UPFC the 
interaction between the real and reactive power flow control system must be reduced. Different control techniques 
for the UPFC system have been proposed. This paper begins with PI decoupling control (PI-D) and it will go on to 
a Generalized predictive controller with neural networks to control the UPFC (NNGPC) [7], [8]. 
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3.1. PI decoupling control (PI-D) 

The principle of this method is to transform the measured magnitudes of the current and the voltage 
of the three phases on the two d-q axes using the park transformation. Next, the values of the active and 
reactive powers are imposed and the reference currents are calculated from these values (the desired powers) 
and the values of the voltages measured, by the two equations: 


3 š $ 
P= 2 (Vsaisa + Vsq isq) (4) 


Q= = (Usa isq — Vsq isa) (5) 


with lra = lä + ipa and irq = isq + ing 


x 2 P*Vsq—-Q"Vsq 

=a a) © 

wx 2 P*Vsq-Q*Vsa 

Re aay aa o 
with A= V24q +V? eq (8) 


where the * superscript defines the reference quantities. The power flow control is then realized by using 
properly designed controllers to force the line currents to follow their references. It is desired that the UPFC 
control system have a fast response with minimal interaction between the real and reactive power together 
with a strong damping of the resonance frequency. However, perfect decoupling is difficult to achieve with a 
PI controller due to the presence of time delays and other non-linearities in the UPFC system. Figure 3 shows 
a step response of decoupling control system. It can be seen that the interaction between the isa and isq appears 
and the system has a slower response. 


PI decoupling control 


P(W) and Q(YAR) 


as 
Time($) 


Figure 3. Power step response of PI with decoupling control system 


3.2. Generalized predictive control (GPC) 
3.2.1. The principle of generalized predictive control (GPC) 

The principle of GPC owes its name to the fact that it takes into account predictions of process’s 
future outputs and sometimes also of the future inputs [9]-[11]. The outputs and the inputs are predicted over 
a finite horizon, then on the synthesis of a control system capable of anticipating setpoint variations. The 
prediction model used is controlled autoregressive integrated moving average (CARIMA), which is an 
extension of the CARMA model, where an integral effect is incorporated to eliminate the permanent 
deviation from the effect of constant disturbances. 


3.2.2. The prediction model and the cost function 

In predictive controllers, several models of the process can be applied, and the implementation of 
the GPC is carried out from the model represented in the form controlled autoregressive integrated moving 
average (CARIMA) Figure 4. He is given by (9). 
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A(ZZ~")y(t) = B (Zut — 1) + x(t) (9) 


x(t) 
Z1B(Z") 


u(t) ve) 


Figure 4. Representation of the CARMA digital model 


With: y(t) output from the system, u(t) command applied to the system, Z7! delay operator, and the 
polynomials are defined by (10). 


ae =1+4Z + + an" (10) 
B(Z-*) = bo + BZ" 4 + bnp? ™ 
x(t)term related to disturbances, usually chosen in the form (11). 

x(t) = Ce eg) (11) 


Hence C(Z~*) = 1 + Z771 +++» + CycZ~"* and e(t) uncorrelated random sequence centered. By 
combining with (9), we obtain the CARMA model: 


A(Z~*)y(t) = B(Z“*)u(t — 1) + C(Z“"*)e(t) (12) 


this model may be inappropriate in several industrial applications in which non-stationary disturbances and in 
this case, the correct model is: 


x(t) = C(Z~*)e(t)/A(Z™) (13) 


A(Z~*)is the differentiation operator which equals: A(Z~*) = 1 — Z~* coupled with (9) gives the controlled 
autoregressive integrated moving average (CARIMA) model as (14). 


A(Z~")y(t) = B(Z*)u(t — 1) + CZ )e(t)/A(Z™) (14) 


For simplicity, we take C(Z~1) = 1 which is the case of white noise, (for the case of C(Z~+) # 1, 
the model becomes (15). 


A(ZZ~")y(t) = B(Z~*)u(t — 1) + e(t)/A(Z*) (15) 


This model constitutes the basic model of the GPC method from which the expression of the control law will 
be derived. To determine the control to be applied to the system at time t, the GPC method minimizes the 
following criterion: 


JO =E n lve +) —we + NP +A n Mult + 7 - 17? (16) 
with Au(t + j) = 0 for j = N, or: N,: is the minimum prediction horizon, N,: is the maximum prediction 
horizon, N,,: is the command horizon, A: is a control signal weighting factor (A > 0), w(t + j): is the 
reference future trajectory known in advance, and Au(t + j): is the control increment. 


To simplify the calculation of the GPC control law, we will set: C(Z~*) = 1, we show that the 
predictor is written: in matrix form: 


J=G.0+F (17) 
where: ¥, U, and F are vectors of dimension NXN with: 

J = [yt + 1), y(t + 2); y(t + NDF (18) 

U = Au), du(t+ 1), n.. ,Au(t + N, — 1]? (19) 
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F = [ft + 1), f(t + 2);...... ft +N (20) 


where: G represents the set of terms involved in the future commands and where f represents the set of 
command terms and past outputs. 


Jo 0 . O 

T 0G 21) 
JN-2 : . 0 
Gn-1 9Jn-2 + Jo 


3.2.3. The predictive control law 

The control law is calculated to minimize the quadratic criterion J (16) and with Au(t + j) = 0 for, 
j 2 N,. The term [y(t + j) — w(t + j)]* represents the weighted sum of the future errors between the future 
outputs and the setpoint signals w(t +j) G =1...... Nz). The term A [Au(t + j — 1]? represents the cost of 
the control effort. J can be written in matrix form as (22). 


J=(V-W)'.(Y —W) +4.U7.U (22) 
The predictor relation is: Y = G.U + F replacing the latter in equation we obtain: 

J=(G.U+F —W)'.(G.U+F —W)+A.U7.U 
J is minimal ire = 0, which gives: 

U = [(G7.G) + AD]"1.G".(F — W) (23) 


where I represent the identity matrix. 

This last equation provides the future increments of control for instants 't' at (t+N,-1), based on the 
information available at instant 't'. Only Au(t) will be applied to the system and the command is then such 
that: w(t) = Au(t) + u(t — 1). The minimization of the criterion will be repeated at each sample for the 
calculation of the new command to be applied to the system Figure 5 [12]-[16]. Or: (yr = i*sa ; u=Vea ; 
(estimated parameters +controller)=minimization of J and y=isa). For the transfer function of our UPFC 
system, the predictive current regulator is synthesized for the adjustment parameters of the Ni, N2, Nu, and A. 
The objective of predictive control is to compute the future control sequence u(t), u (t+1),......... in such a 
way that the future system output y(t+j) is driven close to w(t+j). 


u(t) = (1 —b,b,(b,? +4) u(t — 1) + bybo (by? +A) u(t — 2) +b, (b,? +4) c(t) + 


=A =f 
bi (a, —1)(b,? +a) y(t) -— ba (b? +4) y(t -1) (24) 
estimated parameter y 
parameters estimator 
SE 


u 
na . u UPFC system , 
Controller F = = B a 


Figure 5. Block diagram of the GPC system for UPFC system 


3.3. Neural network generalized predictive control (NNGPC) 

This study seeks to obtain data that will help to address these research gaps for the application of the 
neural generalized predictive control (NGPC) [17]-[21] we follow two steps the first is the identification of 
the system and the second is the design of control. So, in the first step, you develop a neural network model 
system you to control as can you want to see in Figure 6. So, u is our control parameter and it is input to our 
system which will give some output. 
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Neural Network 
Model 


Learning 
Algorithm 


Figure 6. Neural network model predictive control system identification (NNMPCSI) 


We can create this neural network then the error is calculated using the usual method and the apply 
propagation using any of the learning algorithms that we have already chosen. As you can see this is similar 
to our time series neural network so this is the complete model that will be looking into so in the first step, 
we have created the NN model of our system and in the second we will be using this neural network model 
(NNM) to predict the future value. 

The optimization block is also which is the reference signal which we want from the system then the 
optimization block uses some optimization algorithm that optimizes the future system performance. So 
basically, we will have our output from a system which matches nearly with the reference signal Figure 7. In 
this article, the abbreviation will be used to refer to our work. So, these two are the same, one is input to 
NNM, and the second is input to the system to predict. 

The optimization of the controller required the significance of online computation because the 
optimization is performed at each single time step to compute the optimal control. The optimization block 
tries to minimize the performance which is shown in (25). 


I(t) =E en Drt +) Ym lt + HAE n Hj- 1) -— ae j- 2) (25) 
Where: 
Ni, N2, Nu : Horizon over which the tracking error and the control increment are evaluated. 
ú : tentative control signal 
Yr : Desired response 
Ym : neural network model response 
a : determines the effect of control increment on performance 


Data for this study were collected using simulation. Coming back to MATLAB type in the 
command window, write predictor and press enter this will open up the Simulink model Figure 8 [22]-[25]. 
So first let’s look at plant we are using. So, a plant is a set of the equation that represent a dynamic model of 
system (UPFC). 


Controller 


Fr 
Optimization 


UPFC System 


Figure 7. Neural network model predictive control system (NNMPCS) 
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NNG Predictive Controller 


Input 


UPFC System 


Figure 8. NNGPC MATLAB model with system (UPFC) 


This cost function minimizes not only the mean squared error between the reference signal and the 
UPFC system model, but also the weighted squared rate of change of the control input. When this cost 
function is minimized, a control input that meets the constraints is generated that allows the system to track 
the reference trajectory within minor tolerance relying on both configured search parameter (a) and non- 
configured neural network modeling error (e). The following set of J ahead optimal predictions can be 
written by (17) as (26) and (27). 


9 =G.Au(t) + f (26) 
A(t) = bi (b? +4) (c - f) 27) 


Notice that only the first element of u is applied and the procedure is repeated at the next sampling 
time. The algorithm to obtain the control law described in the previous section will be used on the neural 
networks generalized predictive control (NNGPC). Obtaining numerical results for the parameter values 
a, = —0,9231 , b, = 9,6065 1074 , and b, = 2,1617 1077 the horizons being: N, = 1, N, = 2, N, = 1 and 
A = 107°, and to know the function of the control signal for the desired reference with the previous inputs 
and outputs, which are given by a numerical function or in the form of an array: 


u(t) = q u(t — 1) +a, u(t — 2) + æ c(t) + a, y(t) + as y(t — 1) = [a@][u] (28) 
(28) can be expressed as (29). 
u(t) = wy + Wo (29) 


It is known that generally recurrent neural networks (RNN) are considered more suitable for 
modeling dynamical systems and their choice allows us that the training of the network consists in modifying 
the weights and the biases to minimize the quadratic errors in output by using the law of Windrow-Hoff. So, 
with each step of training, the error at the output is calculated as the difference between the required target y 
and the estimated output y of the network. The quantity to be minimized, with each step of training k, is the 
variance of the error at the output of the networks. 


y = wP + wo with wọ =0 (30) 


From which by analogy we obtain: 


w = [a a2 a3 aya | (31) 
u(t — 1) 
u(t — 2) 
P=] c(t) (32) 
y(t) 
yt-1) 


so, P, w and wọ design, respectively, the input vector, the weight, and the bias. Learning takes place 
according to the Widrow-Hoff law. therefore, remember that at each learning step, the output error is 
calculated as the difference between the target sought and the network output. 
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4. SIMULATION RESULTS 

The study was conducted in the form of a survey, with data being gathered via MATLAB/Simulink. 
In the simulation below, in MATLAB neural network and control toolboxes were used and are carried out on 
a Pentium PC under Simulink with model data are chosen equal to the parameters of a laboratory UPFC 
model (Table 1). The NNGP Controller is implemented as C-coded S-functions as shown in Figure 8, and of 
course, with a careful choice of a controller for the inverters. 

Figure 9 shows the step response of the GPC system. It can be seen that the interaction between the 
real force and the interaction with the signals is a bit oscillatory with almost zero prediction error Figure 10. 
Which forces us to apply the recurrent type neural networks to the generalized predictive control for 
improvement Figure 11. So it can be said that the dynamic performance of NNGPC is much better than GPC 
in rising time, set time and overshoot for only in which parameter was tested following a step signal of 
powers (real and reactive). The learning gain of the neural network is fixed at the value of the prediction error 
between the output of the process and the setpoint saturates around a not very large value during a sudden 
variation of the latter, it is equal in absolute value at 0.0258. To verify the effectiveness of the designed 
controllers, different test cases were carried out. 

The test case 1 (robustness test), the inductance of the transmission line is increased by +35% 
compared to the real system. The response of the active and reactive powers to +%35 of XL could not be 
detected, the error message given by “MATLAB command” indicated the saturation of the command to 
infinity during sudden variations of the setpoint signals Figure 12. In test case 2 (stability test), we can check 
the servo response not only in tracking but also in regulation by adding a disturbance. We simulated this time 
by introducing a disturbance of duration 0.02 s and amplitude 10 to once again test the robustness of 
NNGPC. We note that the disturbance does not influence the responses of the powers, on the other hand, its 
action is slightly felt in Figure 13, where the rejection is done gradually and is accompanied by a slight 
overshoot at the external disturbance. 


Table 1. The performance of a Laboratory UPFC model 


Parameters Value Parameters Value 


Ve 220 V Tp 0.80 
V, 220 V Lp 10H 
Vea" 280 V r 0.4Q 
c 2 mF L 10H 


Real power and Reactive power with GPC 


P{W)and Q{YAR) 


o 01 02 03 04 05 06 07 08 09 1 
Time (S) 


Figure 9. Simulation result of step response of system UPFC with GPC 


Number 


Figure 10. Error prediction (GPC) 
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Figure 11. (Real and reactive) powers responses of UPFC with NNGPC 


(Real and Reactive) powers responses at +39% XL (NNGPC) 


o 01 0.2 03 04 05 06 o7 os 0.9 
Time (S) 


Figure 12. Robustness test of NNGPC (at +35% of XL) 


(Real and Reactive) Powers with disturbance at timet = 0.68S (N N GPC) 


2000r T T T T T T T T 


Time disturbance : t =0.02S with amplitude = 10 > 


1000|-- oe ee om 
wv 


= = Į hi a [ T 
1500, 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 


Time(S) 


Figure 13. Stability test of NNGPC (external disturbance) 


ak 


The characteristics of the responses of the UPFC system, Table 2 are deduced graphically from the 
scope block (MATLAB/Simulink). The quality of the adjustment is sometimes specified utilizing certain 
characteristic values of the output of the looped process in slaving. So, the characteristics of the table are 
those of the output of the UPFC system. Thus, for the response time tr (s), we find that it is more or less 
important in the methods using the classical approach than in the methods using neural networks and less 
small in the predictive GPC method using the classical approach than the predictive method using NNGPC 


neural networks. 


Table 2. Characteristics of the different commands 


Command 


t(s) System response time 


t p(s) Disturbance rejection time 


e , Quadratic error 


PI-D 
GPC 
NNGPC 


0.4281 
0.3823 
0.3530 


0.0600 
0.0598 
0.0250 


0.4096 
0.0025, 
0.0016 


Int J Pow Elec & Dri Syst, Vol. 13, No. 3, September 2022: 1448-1458 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 O 1457 


5. CONCLUSION 

This paper proved the quality of robustness and stability characteristics when neural network 
identifier is used with GPC in the UPFC system. It is proved by simulation that NNGPC has satisfied the 
reference pursuit conditions if the appropriate algorithm is used, good offline training is performed, and the 
line search parameter verifies the derived conditions. The simulation results show the quality of the NNGPC 
controller compared to the PI-decoupled controller and the traditional GPC algorithm in terms of meeting the 
required closed-loop performance while stabilizing the system and ensuring robustness. The choice of 
recurrent neural networks with GPC to improve both tracking and disturbance rejection performance of the system. 
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